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The Dependence of Prestellar Core Mass Distributions on the 
Structure of the Parental Cloud 
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ABSTRACT 

The mass distribution of prestellar cores is obtained for clouds with arbitrary 
internal mass distributions using a selection criterion based on the thermal and 
turbulent Jeans mass and applied hierarchically from small to large scales. We 
have checked this methodology comparing our results for a lognormal density 
PDF with the theoretical CMF derived by Hennebelle & Chabrier, namely a 
power-law at large scales and a log-normal cutoff at low scales, but our method 
can be applied to any mass distributions representing a star-forming cloud. This 
methodology enables us to connect the parental cloud structure with the mass 
distribution of the cores and their spatial distribution, providing an efficient tool 
for investigating the physical properties of the molecular clouds that give rise 
to the prestellar core distributions observed. Simulated fBm clouds with the 
Hurst exponent close to the value H — 1/3 give the best agreement with the 
theoretical CMF derived by Hennebelle & Chabrier and Chabrier's system IMF. 
Likewise, the spatial distribution of the cores derived from our methodology show 
a surface density of companions compatible with those observed in Trapezium 
and Ophiucus star-forming regions. This method also allows us to analyze the 
properties of the mass distribution of cores for different realizations. We found 
that the variations in the number of cores formed in different realizations of fBm 
clouds (with the same Hurst exponent) are much larger than the expected root 
Af statistical fluctuations, increasing with H. 
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Introduction 



Is the stellar Initial Mass Function (IMF) universal? This question has been in the lit- 
erature for a long time, and is now extended to the core mass function (CMF) since a close 
relation between the IMF and the CMF has been recognized (Motte et al. 1998; Testi & 
Sargent 1998; Alves et al. 2007; Chabrier & Hennebelle 2010; Michel et al. 2011). Compared 
to the CMF, the mass function of stellar systems seems to be shifted to lower masses by a 
factor that does not depend on the core mass. The currently favored conversion efficiency 
value of the progenitor core mass to the s tellar system is ~ 1/3. However, the origin of this 



conversion efficiency is still contro versial (lAdams fc Fatuzzol Il996l ; iMatzner fc McKed 12000 



Enoch et al. 


2008; 


Dib et al. 


2011) 



ters are observed in the disk of the Galaxy (e.g., iScalol Il998h . but most of these variations 
are consistent with random sampling from a uni versal IMF (jElmegreenl Il997l . Il999i ; Kroupa 
20021 ; iBastian et al.ll2010l ; iParravano et al.l 1201 if ). The non-lineal processes involved in the 
star formation process determine on the one hand the universal form of the IMF and on 
the other hand the range of expected variations of the mass function around this univer- 
sal form. These fluctuation s arise naturally in IMF models based on deterministic chaos 
( jSanchez fc Parrayanolll999n and are also observed in three-dimensional hydrodynamic sim- 
ulations. Recently iGirichidis et al.l ( 120111 ) performed a parameter study of the fragmentation 
properties of collapsing isothermal gas cores with different initial conditions and showed that 
the density profile strongly determines the number of formed stars, the onset of star forma- 
tion, the stellar mass distribution, and the spatial stellar distribution. Furthermore, the 
random setup of the turbulent velocity field in SPH simulation has a major impact in the 
different morphology of the filamentary structure, and consequently on the number of sink 



particles (as shown by IGirichidis et al.l 1201 ll . l2012al jb[). 



The ever increasing resolution of magneto hydrodynamic num erical simulations will pro- 
vide the answer to many of these questi ons (lElmegreen 120111 ). Nevertheless, theoretica l 
I MF models, such as tho s e proposed by iPadoan et al.l (119971 ); iPadoan fc Nordlunl (120021 ) 



or 



Hennebelle fc Chabrierl (120081 120091 ) . provide analytical solutions that help elucidate the 
contribution of the various physical processes involved, but, to obtain these analytic solu- 
tions it is necessary to adopt a series of assumptions that limit their application to specific 
cases. The predictions of these theories have been c ompar ed to the numerical data from 
simulations. In particular, IPadoan fc Nordlundl (12004 ) and IPadoan et al.l (120071) com pared 
the analytical solutions of their theory to numerical simulations. ISchmidt et al.l ( 120101 ) have 
also compared the results from their simulations to the predictions of these theories and have 
shown how the clump mass distribution depends on the turbulence driving mechanism. In 
between these two approaches are the phenomenological models, such as the one presented 
here, that allow one to address some of the questions stated above, in particular, that of 
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sensitivity to the initial conditions. The method consists of a selection criterion based on the 
thermal and turbulent Jeans mass which is applied hierarchically from small to large scales. 



1.1. Aim of the Paper 

The methodology proposed in this work enables a direct connection between the struc- 
ture of molecular clouds and the distributions of generated cores in both mass and space. 
Thus our first aim is to check that the results obtained using our methodology are consistent 
with those obtained using other methods proposed in the literature and that have produced 
reliable results. In particular, we will compare our results for a lognormal density PDF 



with the theoretical CMF derived by iHennebelle fc Chabrierl ( 120081 ). but using two different 



spatial distributions of the cloud mass: a) a ramdon cloud, and b) what we have called 
a "corner" distribution where the voxel mass decreases with the distance to a preselected 
corner. This exercise allows us to evaluate the virtue of the method and how the geometry 
of the cloud defines the dependence of the standard deviation of the lognormal density PDF 
with the smoothing scale R. We have chosen these two very different spatial structures so 
as to make it clear how the analytic formulation of HC08 and the phenomenology presented 
here are connected through the scale dependence of the density PDF. 

Second, we will explore the formation of cores for different parent clouds, but consid- 
ering that the geometry that best describes the spatial structure of the clouds is fractal. 
Observations of close star-forming clouds indicate that the mass distribution in the m can 



be described as having a fractal structure (jFalgarone et al.lll992l ; ISanchez et al.l 120051 1 . The 



analysis will be carried out for fractional Brownian motion (fBm) clouds with a wide range 
of fractal dimensions. Specifically, we will focus on the comparative analysis of the following 
properties: a) Empirical dependence of density PDF on the smoothing scale; b) Mass distri- 
bution of the cores; c) Spatial distribution of the generated cores as measured by the surface 
density of companions, and, d) Cloud core properties averaged over several realizations for 
each H. 

The paper is organized into five sections, this introduction being the first. §2 describes 
the method and defines the main physical variables of the problem and their range of values 
in our simulations. In §3 we check the virtue of our methodology in reproducing the CMF 
derived analytically by Hennebelle & Chabrier (2008), and in §4, we show the application of 
this methodology to fractal clouds generated as fBm clouds with different Hurst exponents 
and compare the results with previous approaches to the same physical systems. Finally, §5 
is devoted to summarizing the main conclusions. 
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2. A Discrete Method for a Hierarchical Collapsing Sequence (HCS Method) 



Following the nomenclature in iMcKee fc Tan! (120031 ) we define a star-forming clump as 



a massive region of molecular gas out of which a star cluster is forming; a core is a region of 
molecular gas that will form a single star (or a multiple-star system such as a binary). The 
resolution at which a distribution of matter is described can be limited by the procedure used 
to generate or measure the distribution, the capacity of storage of information, or simply 
can be chosen to meet a given level of description. In our case, the distribution of matter is 
given in a three-dimensional lattice cube of length L with N 3 0X identical cubic voxels. The 
volume associated with each voxel is l 3 ox , where l vox = L/N vox . The mass of gas contained 
in a voxel centered at coordinates f— l vox x [ix,jy, kz] is denoted as m^j^i where i, j and k 
run from 1 to N vox . 

If we assume that the densest voxel contains a mass m ma ^ and that the physical condi- 
tions in that voxel are such that it is gravitationally unstable □, (i.e. the Jeans length equals 
l vox ), then the thermal Jeans mass M Jyth (oc p -1 ^ 2 ) for a larger cube of d 3 voxels is 

3/2 



\/m cube /d 3 

where m cu i> e = ^ m i,j,k is the mass contained in the volume v whose shape is a cube of 

i,j,k £ v 

side d x l vox . Note that we are implicitly assuming that the temperature and the molecular 
weight are the same in the d 3 voxels in the cube. 

The turbulent Jeans mass can be expressed in terms of the thermal Jeans mass as 

V 3 f R \ 3v d 
M J ,,„ t = M J , ap7 ^(-) =M W -)*, (2) 

where d eq is the length (in l vox units) at which the thermal support and the turbulent 
support are equal, Vo — 1 km s -1 is the turbulent rms velocity at 1 pc scale and C s = 
yj-j^f — 0.22 (T/IOK) 1 ^ 2 (p/2.33) -1 ^ 2 km/s is the sound speed, where 7 g is the adiabatic 
index and /i is the molecular weight. The exponent 77 is the exponent of the Larson's 
(1981) velocity dispersion versus size relation and is related to the 3D power spectrum 
index of the velocity field 77 = (n — 3)/2 w here n = 11/3 for the Kolmogorov case and 4 



for the Burgers case. iKritsuk et al.l (120071 ) estimate n ~ 3.8 — 3.9 (rj ~ 0.4 — 0.45) from 



high resolution hydrodynamic simulations of isothermal supersonic turbulence, in agreement 



1 Larsonl ( 19811 ) noted that "if there is a minimum size of bound condensations produced by supersonic 



compression processes, this may lead to a lower limit of the stellar masses" . 
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with ISchmidt et al.l ( 120091 ) who estimate 77 ~ 0.45 from simul ations of supersonic iso thermal 
turbulence driven by mostly compressive large-scale forcing. IFederrath et al.l (120101 ) showed 
that 77 depends on the nature of the turbulence forcing mechanism; 0.43 for solenoidal forcing 
and 0.47 for compressive forcing. 



Myers fc Fuller! (119921 ) first i ncluded both thermal and nonthermal motions in a model 



of star formation in dense cores. iMcKee fc Tan! (120031 ) focused on the nonthermal part in 
their turbulent core model for massive star formation, but then showed how it is possible 



to smoothly join on to the thermal Jeans mass. Following iHennebelle fc Chabrierl ( 12008 



hereafter HC08), who explicitly included both thermal and nonthermal motions, the Jeans 
mass can be expressed as 

Mj = Mj, th <!!+( — ] \ . (3) 



deq ) 



To apply the Jeans criteria in eq. ([3]) to any cube in the array it is only necessary 
to know friij^ for all i,j, k in the array and the parameters m max and d eq . Note that the 
physical size l vox of the voxels is not needed to determine whether a given cube is Jeans 
unstable. However, as shown below, the parameter d eq depends on the physical conditions 
that determine l vox and on the increase of the velocity dispersion with distance. 

We propose here a procedure to obtain the prestellar core mass distribution which 
is based on a hierarchical collapsing sequence (hereafter the HCS method) in which the 
densest regions collapse first and form the smaller objects. At small scales, thermal support 
dominates and determines the core mass distribution at low masses, whereas, at the largest 
scales turbulence dominates the support and determines the mass distribution at high masses, 
as in the analytical theory of the IMF proposed by HC08. The discrete mass distribution 
rriijji is checked at all scales starting with the smallest, that is at the scale of one voxel, 
i.e. d — 1. By construction only the densest voxel (the one with mass m max ) is marginally 
unstable under the thermal Jeans criterion, however the small turbulent support is enough 
to suppress the collapse at one voxel scale. Then, cubes of side d = 2 are checked to find 
those that fulfill the condition m cu {, e > Mj. In the cubes fulfilling this condition, a core of 
mass m core = Mj is assumed to form giving rise to a stellar system of mass m* = em core . 
The remaining gas m cu b e — e Mj is assumed to become inactive. After this, cubes of side 
d = 3, 4, N vox /2 are considered consecutively. Note that the properties of the velocity 
field are not considered explicitly, but are taken into account implicitly by means of the 
parameters d eq and 77. 

To fix the voxel length l vox we use the assumption that the mass m max contained in the 
densest voxel is gravitationally marginally stable. The radius of a Bonnor-Ebert sphere is 
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R BE ~ 0.486 Rg and the thermal Jeans length is lj )t h = \fn Re, where Rq, the gravitational 
length (McKee & Ostriker 2007), is 

and where njj is the hydrogen nucleus number density. 

Since the density in the densest voxel is uh — 41 ( ^Mq ) / (ipe) cm_3 ; the voxel 
length is 

, , ( m max \ ( 10K n \ 



JM G J V T 2.33, 

where li = 0.15 if lf, ox = ^R BE - Note that this value is close to the value l\ = 0.18 obtained 
if ll ox = |vr (^) 3 . We adopt l x = 0.15 here@ 

Finally, the parameter d eq depends on the turbulent rms velocity which is assumed to 
increase with the size R of the region following the Larson relation (V^,) = V$ x f ^ j , 

where Vo — 1 km s _1 . The value R eq at which thermal and turbulent support are equal can 

1 

be expressed in terms of the sound speed C s as R eg /lpc = (\/3C s /Vo) ''• Therefore, in terms 
of the Mach number at 1 pc scale Mi pc = Vq/C s , the number of voxels d eq for which thermal 
and turbulent support are equal is 

d eq = R eq /l vox = ^Vs/M lpc ^j " /(l vox /lpc). (6) 

We apply the HCS procedure first to mass distributions with a log-normal density 
probability distribution function (PDF) in order to compare our numerical results to the 
analytical CMFs derived by HC08. Later the procedure is applied to fractal clouds with 
density PDFs that are not necessarily log-normal. 



3. Comparison to the HC08 Analytical CMFs Theory 



The analytica l theor y for the IMF developed in HC08 is based on an extension of the 
Press Schechterl (119741 ) statistical formalism applied in cosmology. When applied to the 



2 If we consider that the smallest single brown dwarf that can be formed has a mass of the order of 0.02 
Mq and the progenitor core is about three times this mass, then m max ~ 0.06Af Q . For this value of m ma x 
and the fiducial values T = 10 K and \i = 2.33, eq. © gives l vox = 0.009 pc, a value that roughly agrees 
with the value predicted by Larson's law for this mass 



Larson 1981 ; Kauffmann et al. 2O10f) . 
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mass function of molecular cloud cores, the original Press-Schechter formalism has the prob- 
lem that structures inside structures are not counted. This cloud- in- cloud problem was 
overcome by assuming a conditional probabil ity of finding a collapsed region of mass scale 
M inside a collapsed region of mass scale M' (Ilnutsukall200ll . and references therein). Addi- 
tionally, in the Press-Schechter theory the structures are identified with over-densities in a 
random field of density fluctuations; i.e. a normal distribution in density. Instead, HC08 as- 
sume a log-normal distribution in density, as s uggested by numerical simulations of non-self- 



gravitating supersonic isothermal turbulence (IVazauez-Semadeni 



Passot Vasauez-Semadenilll998l ; l0striker et al 



(IKainulainen et al 



2009 



2001 



1994 



Kritsuk et al. 



Padoan et al. 1997 



20071 ) and observations 



20111 ). Finally, HC08 assume that at any smoothing scale R the 



mass distribution in the cloud is such that the density PDF is always log-normal but with a 
standard deviation <j(R) that decreases with R as 

n*(R) = o>ll-(f) 2n \. (7) 



where Lj is the injection scale and Co is the width of the density distribution at maximum 
resolution R <C Lj. The density PDF at resolution R is then 



Pr(S,R) 



1 



-expy 



2a 2 (R) 



where 5 = ln{p/p) and p is the cloud mean density. For this scale-dependent density PDF 
their theory identifies gravitationally-bound prestellar cores with regions that have a density 
threshold given by the requirement that a fluctuation contains at least one local (thermal 
or turbulent) Jeans mass. As before, the turbulent rms velocity is assumed to correlate 

(jlV 



with size R following the Larson power-law (V^,) = Vq 



The places where the 



average density at scale R is larger than the density threshold contain more than one Jeans 
mass and are expected to form prestellar cores of mass smaller than or equal to the mass 
contained in that region. This is because at smaller scales it may happen that the region is 
not uniform but composed of smaller, denser regions embedded into a more diffuse medium. 
If these denser regions contain one Jeans mass, the end product of the collapse is likely to 
be a cluster of objects whose mass is close to the mass of the smaller /denser regions and not 
to the mass in the volume at scale R. Taking into account the probability of finding these 
unstable sub-structures, HC08 express their core mass function as: 



^HC(^) = 



dM 



oc 



Af tot dlnm (mR 3 ) 



1/2 



1 + (1 - ti)MIR 2t > 
1 + (2r] + l)MlR 2r > 



exp 



[\n(m/R 3 )} 2 a, 



2a 2 



(9) 
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where 



m 



in 
m -. 



(10) 



is the Jeans mass at the average cloud density, R is the radius of the clump in units of 
the Jeans length at the average cloud density, is the characteristic Mach number at the 
Jeans scale, and 77 (~ 0.4 — 0.45) is the exponent of the linewidth-size relation. The ipnci'fh) 
mass distribution in eqs. ([H]- [[(]]) represents the stellar IMF, whereas ^>Hc( m ) represents the 
CMF. 

At low masses (m < 1), the form of ^hc is log-normal. At moderately high masses 



(in > mj) the IMF approaches the power law 



Vhc oc rh-^^ 



rh TiiC 



(m 



> 



1), 



'111 



which gives Thc — 1-3 for the value 77 ~ 0.4 they adopt, in agreement with the iSalpeter 
( 119551 ) value. At very high masses (m 3> 1) the IMF drops off more steeply with mass, 
becoming a log-normal type distribution again. Their results for a non-isothermal equation 
of state are consid erably more complicated, but they are qualitatively consistent with the 
isothermal theory ( jHennebelle fc Chabrierl 120091 ) . 



Note that only the scale-dependent log-normal PDF of the density and the properties of 
turbulence are considered in the HC08 theory. The detailed spatial distribution of matter is 
irrelevant, even when, implicitly, their results refer to the kind of gas distributions associated 
with turbulence. Instead, the HCS method proposed here explicitly takes into account the 
spatial distribution of the matter, and as shown below, the resulting core mass distribution 
is sensitive to this distribution. 



For the moment we do not consider the radiative feedback (Ho 



Gorti k Hollenbach 2002; Krumholz et al. 2007, 2011; Bate 2009 



lenbach fc Tielens 



2012 



Price et all 120091 ) 



1999 



but it is expected that its efficiency also greatly depends on the spatial distribution of the 
matter. 



3.1. Artificial mass distributions with log-normal density PDFs: "corner" and 

"Random" clouds. 

To compare with HC08 we first consider two extreme mass distributions that are useful 
for showing the importance of the spatial structure of the cloud on the resulting mass function 
of collapsing cores. The first is a "random cloud" in which the masses m(i,j,k) and the 
positions k) are uncorrelated. The second is a "corner cloud" in which the locations of 
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the mass voxels m(i,j, k) are ordered in such a way that the density decreases as the sum 
i + j+ k increases. 

Since the lattice is regular and the voxel size l vox is fixed by m max , the mass distribution 
of the voxels (that is the number of voxels with a given mass) is proportional to the density 
PDF, and therefore the standard deviation <jq is the same for the PDF of the masses m(i, j, k) 
and for the PDF of the gas density, and has the same meaning as in eq. ([9]). The PDFs of 
the masses m(i,j, k) in both types of clouds follow the same log-normal function 

N 3 [Mm/rh) - 4? 

dN/dlnm = — v -^exp[- V ' — 2 — , (12) 
cr \/27r 2(T 

where m is the mean mass per voxel and the peak of the distribution occurs at ln(mo) = 
ln{m) — <7q/2. However, the dependence of the standard deviation a(R) on the smoothing 
scale R is, as shown in Fig. [TJ very different. The smooth gradients in the mass distribution 
of the "corner cloud" produce a cr(R) that is close to the HC08 dependence in eq. (JTj). In 
the "random cloud" case, u(R) rapidly drops to zero due to the unphysical discontinuous 
densities that make the average density in any volume containing a relatively small number 
of voxels very close to p. 



3.2. Comparison of HC08 and HCS CMFs for equivalent clouds. 

To compare the core mass function corresponding to a particular mass distribution m^fe 
with ?/>hc we have to determine the set of HC08 parameters (<To, m°, M.*) from the HCS 
input parameters (N mx , m max , T, fi, rj). 

The voxel size l vox (m max ,T, p) is calculated from eq. fl5J), and the mean cloud density 
is n H = m i,j,k/ ( m H L 3 ) with L = N vox x l % 



' vox ■ 



Following HC08 (see also iHennebelle fc Chabrierll2009i ) the Mach number at the Jeans 



scale at the cloud mean density n# is 

Ml ~ (T/lOKf- 1 (r%/10 4 cm- 3 p (13) 

and the Jeans mass at the cloud mean density n# is 

m ° = mx (T/10K) 3/2 (/i/2.33r 3/2 (n^/lOW 3 )" 172 M , (14) 

where the value of m\ depends on the definition of the Jeans mass. HC08 adopt m\ ~ 1, 
but if mj is assumed to be the mass in a sphere of diameter Ij^h = V^Rg then m\ ~ 6.6. If 
instead is assumed to be the mass in a sphere of diameter Ij^jl then m,\ ~ 0.8, close to 







Random 



.1 .2 .3 .4 .5 

R/L 

Fig. 1. — Dependence of the standard deviation on the smoothing scale for two different mass 
distributions and the dependence assumed by HC08 in eq. (JTj) with Lj = L and o"o calculated 
from eqs. ([15] - [16]). The parameters in the three cases are N vox = 2 s , m max = O.O7M , 
T = 10K, \i = 2.33 and rj = 0.4, corresponding to the third case in Table 1. 
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the value adopted in HC08. We assume mi = 1.08, corresponding to the case in which Wj 
is assumed to be the mass in a Bonnor-Ebert sphere. 

The Mach number at the cloud scale Ai = ^-{-r-) v is 

C s K lpc I 

M ~ (1/0.22) (T/10K)" 1/2 (/i/2.33) 1/2 (L/lpc)" , (15) 
and the width of the density distribution o"o at the cloud scale is 

a 2 = ln(l + b 2 M 2 ), (16) 

where b 2 ~ 0.25 (see also Hennebelle & Chabrier 2009). 

For a given set of parameters <To and m max , the masses in the N^ ox voxels in the 
"random" or the "corner" clouds can be assigned following the log-normal distribution 
in eq. f lT2|) . requiring that the mass in the densest voxel is m max . Then the resulting 
core mass function for m^-^ can be compared to V'ho However, what is the appropri- 
ate value for m max to make the comparison? We assume that the appropriate value of 
m max is the value for which the mean cloud density n# is such that the whole cloud is 



close to virial equilibrium; that is, the virial parameter ( )L arson! 11981c iBertoldi fc McKee 



1992|), defined as a VIR = 5a 2 el R/(GM d ) ~ 5(f ) 2 (L/2)/(GM d ) is equal to one. Therefore 



(M d /1M ) ~ 150(L/lpc) 2,?+1 if (Af/lkms^ 1 ) ~ (L/lpc) v . The cloud mass M d = J2 m i,j,k 
depends on the form of the density PDF, on the number of voxels in the lattice and on m max . 
The lattice size L depends on m max through l VO x{jn m ax) in eq. ([5]). 

Table 1 gives the m max values that fulfill the above two conditions when the PDF of 
the mass in voxels is log-normal. Table 1 also gives various derived quantities and the 
corresponding HC08 parameters. Note that the virial parameter is calcu lated omitting the 



effect of the pressure produced by the medium surrounding the cloud. iKainulainen et al. 



(120111 ) estimate that the pressures supporting the clumps against dispersal amount in total 
to about one third of the pressure driving their dispersal. Therefore, cloud masses in Table 
1 exceed the Jeans masses in eq. [3] since M J<th is based on the stability of a Bonnor-Ebert 
sphere. 

For the "corner cloud" and two array sizes, N vox = 128 and 256, Figure [2] shows the 
histograms of the PDF of voxel masses uiij^ an d the core masses obtained with the HCS 
method. Figure [2] also shows the analytical PDF and core mass function from HC08; the 
dashed gray curve is ?/>hc from eq. (Q as function of m = m/mj, whereas the continuous 
gray curve is not normalized to mj. Agreement with the non- normalized HC08 CMF is good 
for all cases in Table 1 . The rapid falloff of the core mass function at low m is due to the fact 
that by construction we have set the most massive voxel to a mass of about 0.06, and that 
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TABLE 1 



Properties of clouds with log-normal density PDFs 



TV 


T 






M cl 


L 


n H 


0o 


M 


M* 


m j 




[K] 


[M e ] 




[M Q ] 


[pc] 


[cm -3 ] 








[M e ] 


64 


10 


0.109 


5.67 


156 


1.05 


5616 


1.42 


5.09 


1.12 


1.44 


128 


10 


0.088 


7.02 


370 


1.69 


3166 


1.53 


6.17 


1.26 


1.92 


256 


10 


0.070 


8.83 


862 


2.69 


1831 


1.64 


7.43 


1.40 


2.52 


256 


8 


0.042 


8.91 


521 


2.02 


2622 


1.64 


7.40 


1.40 


1.51 


256 


12 


0.106 


8.79 


1301 


3.39 


1375 


1.64 


7.44 


1.41 


3.83 




Fig. 2. — Histograms of the PDF of mass in voxels and the resulting core mass function for 
the "corner cloud" case for two array sizes, N vox = 128 (a) and 256 (b). The parameter 
values are those in the second and third lines of Table 1. The continuous gray curve at the 
left of each panel is a log- normal with the appropriate values uq and m max . The continuous 
gray curve at the right of each panel is ipHc( m ) an d the dashed gray curve is the normalized 
HC08 IMF ipK C (m/m°j). 
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sets the minimum mass of a core. At high m the rapid falloff is because the highest mass of a 
core is limited by the mass of the cloud that remains after the formation of the smaller cores. 
Results very similar to those in Fig. [2] are obtained for spherically symmetric distributions 
in which the densest voxel with mass m max is located at the center of the array and the 
mass of the remaining voxels decrease from the center n accord with the log-normal PDF. 
These results indicate that the HCS numerical method captures the main features of the 
HC08 analytical theory. However, the HCS method is not restricted to mass distributions 
with density PDFs following eqs. (J7HS]). 

HC08 express their stellar IMF in terms of which is about 3 M Q for a mean cloud 
density nn ~ 1000 cm -3 . Instead of normalizing masses to we assume that the non- 
normalized function ipnci 171 ) represents CMF and that the mass function of stellar systems 
is shifted to lower masses by a factor e of about 1/3. The evolution of the CMF to the IMF 
has not been definitively established but magnetically-driven outflows (Matzner & McKee 
2000) are expected to produce a mass independent efficiency factor in the range 30 — 50%; 
we assume a value of e ~ 1/3. Note that our procedure for obtaining the stellar mass 
function for a parti c ular g as distrib ution, as well as the theories of the IMF such as those of 
Padoan fc NordlurJ J2OO2I ) (see also lPadoan et alil2007h and HC08, predict the system IMF, 
whereas observations that cannot resolve close binaries determine an effective IMF since 
unresolved binaries are counted as single stars with an effective mass. The determination of 
the individ ual star IMF, in which each star that is a member of a multiple system is counted 
separately (jParravano et al.ll201ll ). is beyond the scope of the present study. The fraction of 
the cloud mass that eventually becomes collapsing cores (~ 0.6), as well as the number of 
these cores, depends weakly on the gas temperature but the mean core mass is proportional 
to (T/10K) 1+1 / , '. 

Contrary to the "corner cloud" or the spherically symmetric distributions, a "random 
cloud" that has the same density PDF does not produce low mass collapsing cores, showing 
that the spatial structure of the cloud is very important not only for the spatial distribution 
of the protostellar objects, but also for the efficiency and mass function of the star forming 
objects. For the PDF of mass in voxels in Fig. (|2]-b), the lowest mass core that collapses in 
the random cloud is 1.3M Q ~ 20m max , whereas for the corner cloud it is O.12M ~ 2m max . 
If the virial parameter of the cloud is increased to olvir = 2.5 (i.e. M c \ ~ 1000M Q and L ~ 5 
pc), then the random cloud produces only a couple of high mass cores (m ~ 100M Q ) but 
the corner cloud still produces a similar number of cores as in the otvm = 1 case. 

It is important to point out that these results show that the dependence of the dispersion 
of the density P DF on the smooth i ng sca le strongly affects the resulting mass distribution of 
collapsing cores. iFederrath et al.l (120101 ) showed that the density PDF in their simulations 
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are roughly consistent with log-normal distributions for both solenoidal and compressive 
forcings, even when the distributions clearly exhibit non-Gaussian higher-order moments. 
However, the dispersion of the density PDF is highly sensitive to the turbulence forcing, and 
therefore they conclude that the theoretical CMF/IMF derived in HC09 is strongly affected 
by the assumed turbulence forcing mechanism. In the following we consider fractal mass 
distributions with roughly log-normal distributions, but with adjustable u(R) functions. 



4. Fractal Clouds 



Frac tal clouds are known to be a good representation of star forming regions (jSanchez et al 



2007al lbl: 1 Efmegreenl |2002| . |2010| ). These kinds of clouds are easy to construct by means of 



recurrence procedures that produce hierarchical self-similar mass distributions. Fractal dis- 
tributio ns are observed over a wide range of scales, from dense cores to giant molecular 
clouds (jBergin fc Tafallal 12007). In particular, numerical simulations of supersonic isother- 



mal turbulence ( Kritsuk et al. 



2007 



Federrath et al.ll2009l ) showed that the density field has 



an approximately fractal structure. 



We focus on fract i onal Brownian motion clouds (hereafter fBm clou ds; IStutzki et al. 



(k x ,ky,k 

ikh 



19981 : lElmegreenl 12002 : iMiville-Deschenes et all 12003 ; ISanchez et all bold ) that have been 
used to represent the internal structu re of molecular clouds. The fB m clouds are generated 
following the procedure described in IMiville-Deschenes et al.l (120031 ) . That is, a field with 
Gaussian distribution intensity Iij^ is obtained by first filling a lattice in wavenumber space 
) with a random phase and Fourier amplitudes proportional to | k \~( H+E / 2 \ where 
ik x 2 + k y 2 + k z 2 y/ 2 , H is the drift (or Hurst) exponent, and E is the dimension of the 
lattice. Subsequently, an inverse fast Fourier transform is applied to generate an intensity 
distribution in real space with a Gaussian intensity distribution. Since intensities must 
be real values, the Fourier amplitudes and phases in the w ave number space have to match 
the appropriate symmetry conditions (IStutzki et al.lll998l ). The exponent H corresponds 
to a power spectrum of the intensity distribution 7 = 2H + E. High resolution numerical 
experiments of s upers onic isothermal turbulence driven by solenoidal forcing, obtained by 
Federrath et al.l (120091 ). have been characterized by a Hurst exponent H = 0.39, correspond- 
ing to a box counting dimen sion Dh — 2.61. Even when th e three dimensional density and 
velocity fields are correlated (iMiville-Deschenes et al.ll2003l ). their relation strongly depends 
on the physical processes involved in the g as dy namics, for example on the energy injection 



mechanism as shown in iFederrath et al.l (120091 ) . HC09 theory takes the exponent of the 
power spectrum of the density and velocity fields into account, but in their CMF analytical 
solutions the two exponents are assumed to be equal. This assumption is used in the previous 
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section to compare our method with HC09 (i.e. through the exponent r\ in eq. ([7])). Note 
that if the two exponents are assumed to be equal, H = 1/3 for a dissipationless cascade of 
energy through an incompressible fluid (i.e. 7 = 11/3 for E — 3). In this section, for the 
application of our method to fBm clouds, we use the quantities r\ and H to parameterize 
separately the velocity and the density fields; the first parameterize the turbulent support, 
and the second determine average properties of the mass distribution in fBm clouds, includ- 
i ng a(R). FBm clo uds have a two-point correlation function of the intensity G/(A) oc X 2H 



( jStutzki et al.lll998l ). so that the mean variation over distance A is Aj (A) oc A 



11 



To generate a t hree-dimen s ional turbulent fractal cloud with a log-normal density dis- 



tribution we follow lElmegreenl (120021 ) : that is, the Gaussian intensity distribution is 
exponentiated to generate a density distribution with a log-normal PDF. The masses mi,j,k 
in the lattice are vn i ~\ t = m max exp\a (7^- 1)1 , where I max is the maximum value of intensi- 
ties Iij,k and a is the contrast factor. Since for a given a the intensity I%j,k is proportional to 
ln(mij,k), the mean variation of ln(m) over distance A is proportional to A . That is, small 
values of H produ ce very rough str uctures, whereas very smooth structures are produced 



with H close to 1. lElmegreenl ( 120021 ) restrict their results to the case H = 1/3. 



By construction the densest voxel has a mass m max . The total mass in the lattice M c i, 
and therefore the mean mass per voxel m, depends on a and H. However, due to the finite 
size of the lattice, two realizations using the same parameter values a and H but different 
randomization of the phases in general do not contain the same total mass M c i. Also, due to 
the finite size of the array, the density PDF departs from the log-normal form, especially at 
low densities and large values of H . To compare the results among fBm clouds we adjust the 
value of a in each realization in order to have the same mass M c i and the same maximum 
voxel mass m max in the lattice. 

Figure ([2]) shows a 2 (R/L) for fBm clouds with various values of the Hurst exponent H. 
Each curve represents the average over 10 realizations and the error bars indicate the 10 and 
90 percentiles. As in Fig. [1] the cloud mass M c i and the value of o"o in the HC08 theory 
(eqlZJ) are the values corresponding to the case in the third line of Table 1. The dashed 
area represents the values of a 2 given by eq.([7|) with an injection scale length in the range 
L/2 < Li< 2L. 

Figure (j3J) shows the mass distribution of the cores formed in 20 different simulations 
with the same parameters H and m max , and the contrast factor a adjusted to always produce 
the same cloud mass. The agreement with iptic{ m ) is much better for the case H = 1/3 than 
for H = 1/6. When H — 1/6 there is a clear excess of high mass cores compared to ^hc(^); 
which can be understood in terms of the dependence a 2 on R. As shown in Fig. [3], except for 
R/L < 0.03 the curve a 2 (R; H = 1/6) is below the HC08 curve a 2 (l-(R/ L) 2r >). When a 2 (R) 
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.1 .2 .3 .4 .5 

R/L 

Fig. 3. — Dependence of the standard deviation on the smoothing scale <J 2 {R) in fBm clouds 
with three different values of the Hurst parameter H. Each curve with error bars corresponds 
to the average value for ten fBm clouds with the same H. The dashed area represents the 
HC08 values of a 2 given by eq. ([7]) with an injection scale length in the range L/2 < Li < 2L. 
The curve with the label HC08 corresponds to eq. (JTj) with Lj = L and a calculated from 
eqs. (EES- HE]). The parameters used are N vox = 2 8 , m max = 0.07 M Q , T = 10K, [i = 2.33 
and rj = 0.4, which correspond to L = 2.7 pc, M c / = 860 M Q , do = 1.64, as in the case of the 
third line of Table 1. 
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is small, the density PDF is a narrow distribution around p and at these densities and sizes 
only massive cores can form. Note that the curve a 2 (R] H = 1/3) intersects the HC08 curve 
at a larger scale and that is why a ffim cloud with H = 1/3 produces in general more low- 
mass cores than that predicted by the HC08 theory. These results highlight the importance 
of the spatial distribution of the gas on the resulting mass distribution of collapsing objects. 

A complex interplay of physical processes determines the final mass of the stars that 
form from a particular configuration of cores in a star-forming region. However, the scaling 
observed between the CMF and the stellar IMF indicate that some average relations between 
these two distributions can be established. For example, the fraction Fh jIOS of individual 
main-sequence stars that will end as core collapse supernovae or the mean mass m* of stars 
in the IMF can be estimated in terms of core to stellar system efficiency (e ~ 1 /3, Matzne r 
& McKee 2000), and the dependence of the binary fraction on the system mass (lLadall2006l ). 



The stellar mean mass can be approximated as m* « (e/ R s )rh core , where R s is the mean 
number of stars per system in the IMF. In other words, R s is defined as the ratio of the total 
number of stars to the t otal number of systems (single star systems + multiple star systems). 
Parravano et all (1201 ll ) estimate that the mean mass of the obje cts in the individual star 
IMF is m* ~ 0.75M Q , so that, assuming that R s ~ 1.3 (lLadall2006l ). the corresponding mean 
core mass is rh core ~ 3M . Note that this estimate of rh core assumes that the core to star 
efficiency does not depend on the mass of the core or on the number of stellar objects formed. 



The fraction of individual main-sequence stars formed with masses over rrih 
be estimated as 

R s ,hN{m core > R S:h m h /e) 



8M & can 



h.ms 



(17) 

R s Af(m core > m bd /e) 

where R s ^ is the mean number of high-mass stars (m > m/J formed in high mass cores 
ijn core > Rg^mh/e). The values quoted in Table 2 correspond to R s ~ 1.3 and R St h — 2. 
Note that eq.(TT71) assumes that the stellar system formed in a core does not disaggregate. 
Additionally, eq.f jT7|) neglects the high-mass primary stars in systems with companions hav- 
ing masses below m^. However, the error introduced is small because most high mass stars 
have a companion of similar mass (Mafz Apellaniz 2008). Eq. f|T7|) also neglects very low- 
mass primaries with brown dwarf companions (m < tum — 0.08 M^) Nonetheless, since th e 
binary fraction of very low-mass stars is small ~ 0.2 (IReid et al.ll2006l ; iBureasser et al.ll2007h. 



the va st majority of systems with m > have main sequence primaries. iParravano et al. 



fl2011f ) estimate that F h:ms ~ (7 - 8) X 10~ 



For seven values of the Hurst exponent H, Table 2 summarizes the average properties 
of the cores formed in sets of 20 simulations of ffim clouds with the same mass M c i but 
a different setup of the random Fourier phases. For the parameter values used for the 
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log m 



Fig. 4. — Mass distribution of the cores formed in a total of 20 different simulations of 
fBm clouds. The solid (dashed) line histogram corresponds to the mass function of fBm 
clouds with H = 1/3 (H = 1/6). The parameter values in all simulations are N vox = 2 8 , 
m max — 0.07 M Q , T = 10 K and \i = 2.33. The contrast factor a in each simulation is 
adjusted to produce a cloud with a mass of 860 M that corresponds to these parameter 
values (see third line in Tablel). As in Fig. (j5J), the sharp rise at core mass of about 
0.1 M is due to the adopted m max value. The continuous curve shows ipnc( m ) f° r these 
parameter values. 
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TABLE 2 





Avera£ 


;e properties of the cores in 


fBm clouds 


H 


A/" a 

* v cores 


f^core 


TP c 

" m,c/g 


TP d 

" h,ms 





28 ±6 


16.28 ±4.90 


0.49 ± 0.06 


0.13 ±0.05 


1/6 


89 ±25 


5.07 ± 1.96 


0.47 ±0.06 


0.022 ±0.02 


1/5 


101 ± 29 


4.47 ± 1.78 


0.47 ±0.07 


0.017 ±0.013 


1/4 


132 ± 38 


3.43 ± 1.43 


0.47 ±0.06 


0.009 ±0.009 


1/3 


201 ±57 


2.42 ±0.87 


0.51 ±0.06 


0.005 ±0.005 


1/2 


405 ± 110 


1.31 ±0.42 


0.57 ±0.06 


0.0006 ±0.0015 


3/4 


797 ± 177 


0.76 ±0.17 


0.67 ±0.05 


0.0000 ±0.0000 



For each value of H, the values quoted correspond to the average over 20 simulations of fBm clouds with 
different random phase but identical parameter values N vox = 2 s , m max = 0.07 M@, T = 10K, = 2.33 and 
rj = 0.4, that correspond to L = 2.7 pc, M c i — 860 Mq, <tq = 1.64, as in the third case in Table 1. The errors 
indicate the standard deviation of the 20 values of each parameter for each value of H 

a Number of cores per cloud. 

& Average mass of cores in a cloud in solar masses. 
c Fraction of the mass of the cloud that collapses in cores. 
^Fraction of stars with masses over 8 M Q estimated with eq. (|17l) . 
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simulations in Table 2, fBm clouds with 1/4 < H < 1/3 produce an average core mass 
distribution that agrees both with the theoretical HC08 CMF and with the expected values 
of fh core pa 3M and -F)j,ms — (7 — 8) x 10~ 3 . Note also that as H increases the number of 
cores increases and their average mass decreases. Figure (J5J) shows the dependence of the 
relative variation of the number of cores AJ\f cores /Af cores as function of H, where AAf cores is 
the standard deviation of M C ores in the 20 simulations. Except for H = the variations in 
the number of cores largely exceed the expected statistical variations. The analysis of 
the dependence of these variations on the considered physical processes and cloud structure 
is out of the scope of the present study. However, we notice here that for a fixed volume of 
simulation and a constant cloud mass the relative variation DeltaN/N is about constant for 
fBm clouds with H values in the range considered in Fig. 5. 



4.1. Spatial distribution of cores 

The spatial distribution of the collapsing cores can be characterized by the surface 
density of companions (SDC) measured on a 2d projection of the positions of the cores by 
sampling all pair of cores over bins of separation AR. In order to compare with the SDC 
in young clusters (Simon 1997) we assume that half of the cores fragment to form two stars 
with 3D separations AR < l vox following a probability distribution p(R) oc AR^ 1 / 2 . Figure 
(|6]) shows the surface density of companions for a single simulation of a fBm cloud with 
H=l/3 and the same parameter values used in Fig. (j3J) and the cas e H=l /3 in Table 2. 



The labeled power-laws in Fig. (JBJ) are the SDCs reported by ISimonl (119971 ) for the Orion 
Trapezium star formation region and for the Ophiuchus star formation region. The SDC for 
this particular simulation is in between the observed SDC's for these two regions. The fall 
of the SDC at large radii is due to edge effects when AR is of the order of size L of the 
simulation. For other simulations with the same parameter values the SDCs are similar, but 
when the parameter H is increased the SDC curves shift upwards. 



5. Conclusions 

We have proposed a procedure to obtain the prestellar core mass distribution that results 
from the collapse of prestellar cores in clumps by assuming that the densest regions collapse 
first and form the smaller objects. At small scales, thermal support dominates and deter- 
mines the mass distribution of cores at low masses, whereas at the largest scales turbulence 
dominates the support and determines themass distribution at high masses. The numerical 
method proposed here make use of a small number of parameters, namely N vox , T, fi, r], to- 
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< 1/3 1/2 3/4 1 

en 

H 

Fig. 5. — Relative variation of the number of cores (AJ\f cores /Af cores ) in 20 simulations as 
function of the Hurst parameter H. The open circles are the corresponding values 1/VJ7 
expected if they were statistical variations of the number of cores. 
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Fig. 6. — Surface density of companions S(Ai2) vs the separation in AU for a single sim- 
ulation of fBm clouds with H=l/3. The small dots connected by dashed lines correspond 
to E(Ai2) for the cores. The large dots connected by solid lines correspond to the surface 
density of companions systems (AR > l vox ) when half of the cores are assumed to frag- 
ment into two units. The open circles connected by light lines correspond to the surface 
density of companion stars assuming that the binary separation AR is less than l vox with 
a probability distribution p(R) oc AR^ 1 ^ 2 . The two segment power-law gray lines represent 
the surf ace density o f companions for the Orion Trapezium and Ophiuchus star formation 
regions (ISimonlll997l ). 
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gether with the cloud properties (i.e. M c \ and H) and the assumptions that the mass in the 
densest voxel m max is equal to the Jeans mass and that the cloud as a whole is marginally 
stable. When the proposed method is applied to a mass distribution whose density PDF is 
a log-normal at all smoothing scales R and its standard deviation cr(R) is given by eq. 
the average mass distribution agrees with the CMF predicted by the analytical theory of the 
IMF proposed by HC08. The HCS method can be seen as a numerical version of the HC08 
theory, and there is univocal correspondence between the parameters in both models. 

Both the Padoan-Nordlund IMF and the Hennebelle-Chabrier IMF apply to particular 
star-forming cloud conditions. In order to determine an average IMF that can be compared 
with observations of stars from different clouds, it is necessary to average their theoretical 
IMFs for a distribution of cloud temperatures, densities and Mach numbers. The core mass 
distribution from our method is even more dependent on cloud property since, as we have 
shown, the masses of the resulting cores also depend on the particular distribution of mass 
within the cloud. Large variations in the resulting core mass distribution are observed in 
ffim clouds with the same mass M c \ and Hurst exponent H, but a different setup of the 
random Fourier phases. As shown in Table 2 and Fig. (151), the number of cores in a set of 
20 simulations display variations that largely exceed the expected \0\f statistical variations. 
Due to its simplicity the HCS method is computationally efficient at obtaining the mass and 
position of the cores that collapse in an arbitrary distribution of gas. Therefore the HCS 
method is well suited to analyzing the effects produced by changes in the physics over a large 
number of initial conditions. 

We have applied the HCS method to lattices with a number of cells up to 2 8x3 , which 
represent clumps of mass ~ 10 3 M and size ~ 3 pc, but larger lattices can be processed. 
There is no restriction in the way the mass in the voxels is assigned, but we have focused on 
fBm clouds that have been used as analogs of real interstellar clouds. We confirm that fBm 
clouds with H ~ 1/3, corresponding to 7 = 11/3 (Elmegreen 2002), give better agreement 
with the theoretical CMF derived by Hennebelle and Chabrier and the observed IMF. We 
have also shown that the spatial distribution of the cores for fBm clouds with H = 1/3 has 
a surface density of companions that resembles that of young stellar clusters (Simon 1997). 
Since the HCS method provides the sequence and location of newly formed stars, the method 
can be easily modified to consider radiative feedback effects. 

We wish to thank David Hollenbach and Christopher McKee for their valuable com- 
ments. AP was supported as guest researcher for part of this research by the IESA-CSIC. 
N.S. acknowledges financial support from the Ministerio de Economi'a y Competitividad of 
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